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ABSTRACT 

Context. The properties of the dust grains (e.g., temperature and mass) can be derived from fitting far-IR SEDs (>100 ;um). Only with SPIRE on 
Herschel has it been possible to get high spatial resolution at 200 to 500 yum that is beyond the peak (~ 160 pm) of dust emission in most galaxies. 
Aims. We investigate the differences in the fitted dust temperatures and masses determined using only <200 yum data and then also including 
>200 ytim data (new SPIRE observations) to determine how important having >200 yum data is for deriving these dust properties. 
Methods. We fit the 100 to 350 yum observations of the Large Magellanic Cloud (LMC) point-by-point with a model that consists of a single 
temperature and fixed emissivity law. The data used are existing observations at 100 and 160 yum (from IRAS and Spitzer) and new SPIRE 
observations of 1/4 of the LMC observed for the HERITAGE Key Project as part of the Herschel Science Demonstration phase. 
Results. The dust temperatures and masses computed using only 100 and 160 yum data can differ by up to 10% and 36%, respectively, from 
those that also include the SPIRE 250 & 350 yum data. We find that an emissivity law proportional to /I"' '' minimizes the 100-350 yum fractional 
residuals. We find that the emission at 500 yum is ~10% higher than expected from extrapolating the fits made at shorter wavelengths. We find the 
fractional 500 yum excess is weakly anti-correlated with MIPS 24 yum flux and the total gas surface density. This argues against a flux calibration 
error as the origin of the 500 yum excess. Our results do not allow us to distinguish between a systematic variation in the wavelength dependent 
emissivity law or a population of very cold dust only detectable atA> 500 ^m for the origin of the 500 yum excess. 
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1. Introduction 

Among nearby galaxies, the Large Magellanic Cloud (LMC) and 
Small Magellanic Cloud (SMC) represent unique astrophysical 
laboratories for interstellar medium (ISM) s tudies. Both Cl ouds 
are relatively nearby, the L MC at -50 kpc (ISchaefeiil2008l) and 
the SMC at -60 kpc (Hil ditchetal.ll2005h . and provide ISM 
measurements that are relatively unconfused by line-of-sight un- 
certainties when compared to the Milky Way. The two Clouds 
span an interesting metallic! ty range with the LMC at -1/2 Zq 
jRussell & Dopitalll992,) being above the threshold of 1/3-1/4 
Zq where the properties of the ISM change as traced by the 
rapid reduction in the PAH dust mass fractions and possible 
dust-to-gas ra tios (Drai neet al. 2007) and the SMC at ~1/5Zq 
dRussell & Do pita 199^ below this threshold. Finally, the dust 
in the LMC and SMC shows stron g variations in its ultraviolet 
characteristics (Gordo n et al.ll2003h . 

The HERschel Inventory of The Agents of Galaxy Evolution 
(HERITAGE) in the Magellanic Clouds Herschel Key Project 
will map both Clouds using the PACS/SPIRE Parallel observing 
mode providing observations at 100, 160, 250, 350, and 500 yum 
(Meixner, M. et. al. 2010). The HERITAGE wavelength cover- 
age (100-500 yum) and spatial resolution (-10 pc at 500 jum) is 
well suited to measuring the spatial variations of dust tempera- 
tures and masses. The infrared dust emission in most galaxies 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



peaks between 100-200 yum dPale et al.ll2005h and observations 
>200 yum are important for accurate dust temperature and masses 
(iWillmer et al . 2009). Ground-based submilimeter observations 
do provide the needed >200yL(m observations, but they have been 
seen to be in excess of that expected from extrapol ating fits to the 
<200 jum data for sub-solar metallicity galaxies dGalhano et al.l 
2005). This excess could be due to very cold dust that only emits 
at submilimeter wavelengths o r variations in the wavelength de- 
pendent dust emissivity law dReach et al.l 119951 : iParadis et al.l 
2009a). As part of the Science Demonstration Program (SDP), 
two HERITAGE AORs centered on the LMC were executed. 
These observations are used in this paper to explore the impact 
SPIRE observations have on the measurement of dust temper- 
atures and masses including the behavior of any submilimeter 
excess. 



2. Data 

The observation a nd data reduction f o r the HERITAGE SDP 
data are given in 'Meixner. M. et. al] d201CA. For this paper, 
we use high quality IRAS 100 /im and MIPS 160 //m ob- 
servations instead of the PACS observations which display 
large residual instrumental signatures (expected to be elim- 
inated with the full HERITAGE dataset). We extracted the 
HERITAGE SDP region from the existing IRAS/IRIS 100 /urn 
(Miville-Deschenes & Lagache 2005 ) and MIPS 160 /urn 
dMeixner et alj 120061: [Bernard et"aiT r2008) mosaics. We have 
used custom convolution kernels created using the technique of 
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Fig. 1. The best fit (assuming an emissivity law withyS - 1 .5) dust temperature and mass images are shown along with the fractional 
residual images at all 5 bands. The resolution of all the images is 43. Each of the fractional residual images shows the difference 
in flux between the measured value and best fit model divided by the best fit model. They are separately linearly scaled (black to 
white) to emphasize the structure and the scale ranges are given at the top of each of the images. The dust temperature image is 
linearly scaled between 18 (blue) and 28 K (yellow). The dust mass image has a sqrt scaling between 10^ (black) and 10'' (yellow) 
Mo/sr. The vertical white and black streaks in the SPIRE fractional residual images are caused by residual instrumental signatures. 



iGordon et al.l (l2008l) to convolve the images to a common reso- 
lution of 4^3 of the IRAS 100 /vm data. We also created a 2nd 
set of images (excluding the IRAS 100 /im data) at the common 
resolution of the SPIRE 500 fim and MIPS 160 fim data of -38". 
Emission from Milky Way (MW) foreground cirrus clouds 
contributes to the far-IR emission seen in the LMC. We use the 
HI column density map created by integrating the MW veloc- 
ities in the full HI cube fStavel ev-Smith et al.l *2003) to correct 
all the images for the MW infrared cirrus emission. The HI col- 
umn densities were transformed to IR surf ace brightness e s usin g 
the model of the MW emission used by iBernard et al.l (l2008h . 
Finally, any residual emission was removed by fitting a gradient 
across the SDP region using the regions in the strip beyond the 
IR edge of the LMC. 



3. Results 

For each point in the image, we determined the dust temperature 
by fitting the observed far-IR SED to a modified black body of 
the form 



F, ex A-PB,{Tdus,) 



(1) 



The dust mass is computed from the measured 160 yum flux 
{Fiho), at each point, using 
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(2) 



where a = 0.1 /im is the grain radius, the grains are assumed 
to be spherical silicate grains with a density p = 3 g cm"^, 
d = 50 kpc is the LMC distance, and ee„,(160) = 5.5 x 10""^ 
( Laor & Draine 199 3). This method is fairly standa rd and while 
other more sophisticated fi tting methods exist (iDraine et al.l 
I2007L iGallianoetaLl 120081: iParadis et alJ l2009bl) . this simple 
model allows us to probe the effects of adding >200 ^im data 
to the fits with fewest assumptions. We restrict our fits to us- 
ing only data >100 /^m as observations at shorter wavelengths 
likely include non-equilibrium dust grain emission (transient 
heating). The data points used in the fits were weighted by the 
uncertainties (i.e., 1/cr^). The main uncertainties on the mea- 
surements are the calibration and background noise uncertainties 
and we sum them in quadrature. The calibration errors are as- 
sumed to be approximately the same at all bands a t around 15% 
dStansberrv et alj|2007t [Swinvard. B. et. alj|2010l) . Data within 
Icr of the background are not used in the fits. 



3.1. Dust temperatures and masses 

The best fit dust temperature and mass values were determined 
by for fits using only the pre-Herschel data (IRAS 100 //m 
and MIPS 160 fim) and fits including the Herschel SPIRE data 
(IRAS 100 yum, MIPS 160 yum, and SPIRE 250/350 yum). Given 
the inclusion of the IRAS 100 yum, we used the 4'3 resolution im- 
ages. The SPIRE 500 yum data are not include d in these fits as it is 
usually systematically high (see Sec. l3.2l and lMeixner. M. et. al.l 
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Fig. 2. The histograms of the fractional residuals at different 
wavelengths are shown for/3 =1.0, 1.5, and 2.0. The fractional 
residual is the difference in flux at each wavelength between the 
measured and best fit model divided by the best fit model. The 
jS = 1.5 and 2.0 histograms are offset by 3000 and 6000 pixels, 
respectively. The dashed vertical line indicates zero fractional 
residual. The strongly peaked 350 ;um histogram is simply a re- 
sult of the relative weighting of different wavelengths in the fit. 
For example, using equal weights produces more equally peaked 
histograms between the different wavelengths. 



( I2OIOI) ) and including it in the fits only causes the residuals at 
the other wavelengths to increase without significantly improv- 
ing the fit. The value of/3 used in the fits was set to 1, 1.5, or 2 as 
this range encompasses realistic dust grains (amorp hous to crys- 
talline grains) and is what has been used in the past ( Dunne et al.l 
I2OOOI) . The dust temperature and mass maps and fractional resid- 
ual images for the y6 = 1 .5 case are shown in Fig. [T] 

The differences between the best fit dust temperatures and 
masses (with and without the SPIRE data) depends on the value 
of jS used. ForyS = 2, the with SPIRE to without SPIRE tem- 
perature ratio is 0.97 + 0.06 and mass ratio is 1.19 + 0.31. For 
/3 - 1.5, the with/without SPIRE temperature ratio is 1 .02 + 0.07 
and mass ratio is 0.96 + 0.25. For /J - 1, the with/without tem- 
perature ratio is 1.08 + 0.08 and mass ratio is 0.77 + 0.20. Thus, 
the inclusion of > 200 fim data in the fits can change the derived 
dust temperature by up to 8% and mass by up to 23% depending 
on the assumed value of /3. 

Prior to the Herschel observations, it was not possible to con- 
strain the best value of/? given that there were only two infrared 
maps of the LMC with A > 100 jjm. With the Herschel obser- 
vations, the behavior of the residuals as a function of /? can be 
used to determine the optimal /3 value. Histograms of the frac- 
tional residuals at different/? values are shown in Fig.|2] A value 



of/? = 1 .5 clearly minimizes the fractional residuals at all wave- 
lengths with most of the pixels having residuals of less 10% at 
all wavelengths except 500 //m. This result implies that either 
the characteristics of the dust grains are intermediate between 
the two extremes or that a more complex dust emission model 
includin g a distribution of d ust temperatures an d grain sizes is 
needed ( Draine et al.l 120071 : iParadis et al.l l2009b). Assuming a 
/3 = 2.0 for the pie-Herschel fits (a common assumption) and us- 
ing the best fit jS = 1.5 for the fits including the SPIRE data, we 
find the with/without temperature ratio is 1.12 ± 0.07 and mass 
ratio is 0.64 + 0.16. This decrease in du st masses reduces the 
magnitude of the "FIR excess" found by Bernard et alj (l2008h 
for the LMC. iRoman-Duval. J. et. al.l (1201 Ol) explore this issue 
in detail for two specific LMC molecular clouds. 

3.2. 500 fj.m excess 

In the previous section, we have not included the 500 jum ob- 
servations in the analysis as it was see n not to improve the 
quality of the fits and previous studies ( Galliano et alj 120051 : 
[Galametz et a l. 2009) have observed submm fluxes in excess of 
that expected from fits to the far-IR fluxes. At 4'3 resolution, the 
average fractional 500 //m fit residual is 0.25, 0. 10, and -0.05 for 
/3 values of 2, 1.5, and 1 (Fig.|2]i. As a/3 - 1.5 is strongly favored 
as it minimizes the residuals at all other wavelengths, we find a 
500 jum excess of approximately 10%. We find the same level of 
500 fiui excess for fits done at both the 4'3 and 38" resolutions. 

There are four possible origins of the 500 /urn excess: 1) sys- 
tematics due to our assumptions on our fitting, 2) a flux cali- 
bration err or, 3) variations in the wavelength dep endent emis- 
sivity law dReach et alJI 19951: lAgladze etal]|1996h, and 4) very 
cold dust that mostly emits at > 500 /im dFinkbeiner et al.l[l999l : 
iGalliano et al.ll20()5l) . Whatever the the physical process respon- 
sible for the 500 jum excess, the HERITAGE SDP SPIRE data 
of the LMC allow us to probe the origin of the 500 fj.m excess 
at high spatial resolution in an external galaxy for the first time. 
We tested the systematics of our fitting algorithm and searched 
for correlations of the 500 fim excess with different tracers of 
the ISM conditions (dust temperature, dust mass, HI mass, and 
MIPS 24 /urn flux) in an attempt to determine the origin of the 
500 //m excess. The two strongest correlations are seen for MIPS 
24 yum flux (probing the ISM conditions for small grains) and the 
total gas mass (probing the ISM density) are are shown in Fig. [3] 

To test 1), we performed Monte Carlo simulations where the 
observations were simulated both with and without an excess at 
500 jum and with different /? values. These simulations were fit 
with varying /? laws and realistic uncertainties. A 500 yum excess 
was found in the simulations if it was part of the simulation or 
if the fitting /? was smaller than the simulation /3. Given that we 
empirically determine /? from the < 500 //m data, our conclusion 
is that the excess we find is not a result of our fitting method. 

For 2), it is possible that there is a systematic 500 fim flux 
calibration error on the order of 10%. The official maximal 
possibl e flux calibration error for SPIRE is 15% at any wave- 
length (iGriffin. M. et. al.l[2010l) . Given that we are including the 
SPIRE 250 and 350 yum measurements in our fitting, the 500 //m 
flux calibration error would have to be relative to the other two 
SPIRE bands and so is likely smaller than 15%. In addition, we 
would expect to see no correlation between the excess and ISM 
condition tracers, yet we see weak correlations (Fig.|3]l. 

For 3), a wavelength dependent increase in the dust emissiv- 
ity law at 500 /im on the order of 10% is possible (Paradis et al] 
2009a). This variation may be attributed to the dust grains amor- 
phous/crystalline nature, size distribution, temperature, and ma- 
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Fig. 3. The 500 //m excess is plotted versus MIPS 24 //m flux and total gas mass surface density using the 38" resolution images. 
The total gas mass is deriv ed from HI and CO observations (Kim et al. 2003; Fukui et al. 2008) assuming Xqq = 7 x 10^" cm"^ 
(K km s ')"' dFukui et al.l l2008) and an overall 36% increase in the gas mass to account for the associated He gas. The density of 
points is coded from low to high as purple -red-green-black. The solid and dashed black lines give the zero and 10% excess values, 
respectively. The solid points give the values of the mode in equally spaced logarithmic bins. On average, the residuals are markedly 
negative at low MIPS 24 /im fluxes indicating that below these fluxes, the residuals are affected by background subtraction or fitting 
errors. 



terial (lHenningetal.1119951; lMenvetal.ll2007h . For example, if 
the 500 iiva. excess is due to small grains having a different yS 
than large grains, we would exp ect the excess to b e correlated 
with the MIPS 24 //m emission (Rea ch etaLlll995h . Yet the ex- 
cess is weakly anti-correlated with MIPS 24 /zm flux (Fig.|3]l. 

For 4), very cold dust that only emits at >500 ytva is phys- 
ically possible. The very coldest dust would necessarily be the 
dust that is best shielded and, thus, we would then expect the 
excess to be strongest in the highest density regions and regions 
with the lowest radiation fields. Figure [3] gives a conflicting an- 
swer as we see the largest excesses in the faintest 24 //m regions 
(as expected) and least dense regions (not as expected). 



4. Conclusions 

We investigate the importance of >200 //m data in determining 
dust temperatures and masses using new Herschel SPIRE ob- 
servations of the LMC (taken for the HERITAGE Key Project 
as part of the Herschel Science Demonstration phase) combined 
with existing IRAS 100 ptm. and Spitzer MIPS 160 fim images. 
We fit the observations with a model consisting of dust emitting 
as a single temperature blackbody modified with an emissivity 
law proportional to A^^. For fixed values of/3, fits using only the 
100-160 jjm data give dust temperatures and masses that are on 
average up to 8% and 23% different from fits using the sameyS 
and the 100-350 ;um data. The new SPIRE observations allowed 
us to determine thatyS =1.5 minimizes the residuals from 100 to 
350 /jm. Using ayS = 2.0 for the 100-160 //m and a/? = 1.5 for 
the 100-350 fim fits results in an increase of 10% for the dust 
temperature and a decrease in the dust mass by 36%. 

On average, there is a fractional excess at 500 jim of ~10%. 
The origin of the fractional excess is unlikely to be due to our 
fitting algorithm or a flux calibration error, but it could be due to 
either very cold dust that emits only >500 yum or a variation in 
the wavelength dependent change in the dust emissivity. Planned 
HERITAGE observations of the LMC and SMC will allow for a 
more detailed investigation of including > 200 fim data (mainly 



the 500 /im excess) due to better quality PACS and SPIRE im- 
ages (optimized observations and cross-scans). 
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